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ABSTRACT 

Context. Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray 
detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power 
pulsar (£ = 3.1 x 10 36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise 
typical of very young pulsars. In the very-high-energy domain (E > 100 GeV), H.E.S.S. has reported the detection of the extended 
source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. 

Aims. We search for gamma and X-ray pulsations from this pulsar, characterize the neutron star emission and explore the environment 
of PSRJ1357-6429. 

Methods. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1 .4 GHz, we phase-fold more 
than two years of gamma-ray data acquired by the Large Area Telescope aboard Fermi as well as those collected with XMM-Newton, 
and perform gamma-ray spectral modeling. 

Results. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one 
broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 
1.5 ± 0.3 stat ± 0.3 sysl with an exponential cut-off at 0.8 ± 0.3 stat ± 0.3 syst GeV and an integral photon flux above 100 MeV of (6.5 ± 
1.6 st at ± 2.3 syst ) xl0~ 8 cirT 2 s -1 . The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin 
(2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. 
Conclusions. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of 

L y = (2.13 ± 0.25 st!it ± 0.83 syst ) x 10 34 erg s _1 , consistent with an L y <x relationship. The Fermi non-detection of the pulsar 
wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV 
emission. 

Key words, pulsar, pulsar wind nebula, PSR J 1357-6429 



1. Introduction 

Since the launch of the Fermi Gamma-ray Space Telescope, the 
number of known gamma-ray pulsars has increased dramati- 



cally. T he first Fermi Larg e Area Telescope (LAT) catalog of 
pulsars dAbdo et alll2010al) also marks the birth of new classes 
of gamma-ray pulsars, such as millis econd dAbdo et all l2 009a) 
and radio-faint gamma-ray pulsars dAbdo et all l2009b). The 
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rest are young pulsars, discovered via an efficient collaboration 
among gamma-ray, radio and X-ray astronomers. Interestingly, 
many of the detected pulsars are found to be powering pulsar 
wind nebulae (PWN), and some are associated with TeV sources. 

PSR J 1357 -6429 is a very young, energetic pulsar dis- 
covered during the Parkes multib e am su rvey of the Galactic 
plane, as reported in Camil o et al.l d2004f) . From its dispersion 
measure of (127.2 + 0.5) cm' ] pc, the NE200 1 mod el of the 
Galactic electron distribution dCordes & Lazi assigns 
PSR J1357-6429 a distance of d = 2.4 ± 0.6 kpc. With a spin 
period of 166 ms and a period derivative of 3.6 x 10~ 13 s s _1 , the 
characteristic age is r = 7300 yr. 

No EGRET source is coincident with PSR J1357-6429. 
Recently, Very Large Telescope data were used to search fo r op- 
tical emission from PSR J1357-6429 (Mig nani et ail 1201 lh . No 
counterpart was found, and an upper limit of 27 was de- 
termin ed. Using X-ray Chandra data collected in 2005. IZavlinl 
(2007) reported a hint of pulsed emission using an accurate, 
contemporaneous radio rotation ephemeris. While of low sig- 
nificance (2-3cr), it is nevertheless mor e compelling than th e 
non-detection of pulsations reported by Esposi tcTet alJ (120071) . 
who folded the same data with an older ephemeris as a starting 
point to scan for a rotation period for which X-ray pulsations 
might appear. This pulsar's timing noise imposes a large num- 
ber of period trials, which in turn diminishes the significance of 
what was, at best, a weak signal. In this article, we confirm the 
detection using new XMM-Newton data from 2009. 

PSR J 1 357-6429 was highlighted bv lAbdo et al] (l2010al) be- 
cause of the absence of gamma-ray pulsati ons, and might hav e 
become a candidate "gamma-faint" pulsar (Rom ani et al.Ll201 lh . 
Nearby, with a high spin-down power E = 3.1 X 10 36 erg s _1 , it 

has VE/d 2 that is 1% of Vela's, a threshold above which the 
LAT detects most pulsars. Here we show that the previous non- 
detection was due not to beaming or emission physics, but rather 
because the timing noise typical of very young pulsars made it 
difficult to construct an accurate rotation ephemeris after a large 
glitch early in the Fermi epoch. Once enough post-glitch radio 
observations had been obtained to build a phase-connected tim- 
ing model, the gam ma-ray pulsations beca me immediately ap- 
parent. We note that lPellizzoni et a l. (2009) reported a weak de- 
tection with AGILE. 

Very often, young and energetic pulsars are embedded 
in compact nebulae powered by relativistic pulsar winds. 
Short X-ray observations were performed in 2005 with XMM- 
Newto n and Chandra to search fo r X-ra y emission from a 
PWN (Esposit o et aTl l2007t IZavlinl |2007l) . Using the XMM- 
Newton data, lEsposito et al.l (120071) found marginal evidence of 
diffuse emission in the 2—4 keV energy band consisting of a faint 
elongated structure starting from PSR J1357-6429. They also 
set a 3<x upper limit on the luminosity of ~ 3 x 10 31 erg s _1 in the 
2-10 keV range on putati ve diffu s e emi ssion using the Chandra 
data. On the other hand, IZavlinl d2007l) reported the detection 
of tail-like emission associated with the pulsar with a 0.5-10 
keV luminosity of ~ 2.5 x 10 31 erg s _1 . In the very high energy 
domain (E > 100 GeV), the H.E.S.S. data show significantly 
extended emission (intrinsic Gau ssian width of 0.20° + .02°) 
lying close to PSR J1357-6429 dAbramowski et all 1201 ll) . In 
this article, a GeV upper limit is derived on the emission from 
the possible PWN observed in the TeV band, providing new 
constraints on the leptonic population responsible for the very 
high energy emission. 



2. LAT description and observations 

The LAT is a gamma-ray telescope that detects photons by con- 
version into electron-positron pairs in the energy range between 
20 Me V to more than 300 GeV, as described by lAtwood et al.l 
d2009l) . It is made of a high-resolution converter/tracker (for 
direction measurement of the incident gamma rays), a CsI(Tl) 
crystal calorimeter (for energy measurement) and an anti- 
coincidence detector to identify the background of charged par- 
ticles. Compared to EGRET, the LAT has a larger effective area 
(~ 8000 cm 2 on-axis above 1 GeV), a broader field of view (~ 
2.4 sr) and superior angular resolution (~ 0.6° 68% containment 
at 1 GeV for events converting in the fr ont section o f the trac ker). 
The on-orbit calibration is described in lAbdoetalJ d2009ch. 

The analysis used 29 months of data collected starting 
August 4, 2008, and extending until January 15, 2011. We se- 
lected events with energies greater than 0. 1 GeV but excluded 
those with zenith angles larger than 100°, to minimize con- 
taminatio n from secondary gamma-rays from the Earth's atmo- 
sphere dAbdo et all l2009dl). We se l ected events in the Source 
class, as described in lAtwood et a l. (2009). Source events opti- 
mize the trade-off between gamma-ray detection efficiency and 
residual charged-particle contamination in a manner well-suited 
to localized, persistent sources, as distinct from the Transient or 
Diffuse classes. We used the most recent analysis version, called 
"Pass 7", documented at the Fermi Science Support CenterQ. 

3. Radio timing observations 

As part of the space- an d ground-based pulsar timing 
campaign supporting Fermi dSmith et all 120081) . the pulsar 
PSR J 1357-6 429 is being obs e rved a t the 64-m Parkes ra- 
dio telescope ( Weltevre de et all 120101) . A total of 74 obser- 
vations spanning MJDs 54220 to 55575 were made in a 
band centred at 1369 M Hz. In this band, th e pulsar has a 
flux density of 440 /Jy dCamilo et all 120041) . Observations 
were ma de using the central beam o f the 20 cm Multibeam 
Receiver dStavelev-Smith et all [T996) and recorded using one 
of the digital filterbanks (DFBs) built by the Australia Telescope 
National Facility. The DFBs employed real time pulse phase 
folding to average data into 30 second sub-integrations using 
256 MHz of total bandwidth. After calibration, pulse profiles 
for ~ 5 minute observations were cross correlated with a tem- 
plate to obta i n time s-of-arrival. The TEMP02 timing package 
dHobbs et all 12006) was then used to build the timing solution, 
taking into account a significant glitch at MJD 54786 and strong 
timing noise prevalent in young pulsars. The glitch was mod- 
eled by including jumps in spin phase, spin frequency, and spin 
frequency derivative at the glitch epoch. Higher order spin fre- 
quency derivatives and harmonically related sinusoids (Hobbs et 
al 2004) were used to model the timing noise. The post-fit rms 
is 1 .4 ms, or 0.8 % of a pulsar rotation. This timing solution will 
be made available through the FSSC 1 . 

4. X-ray observations and data reduction 

XMM-Newton observed PSR J1357-6429 on October 14-15, 
2009 with the EPIC-pn instrument operated in Small Window 
mode and two EPIC-MOS detectors in Full Window mode (for 
55.5 and 78.1 ks effective exposures, respectively). The data 
were reprocessed and reduced with the XMM-Newton Science 
Analysis Software (SAS v. 11.0). 

1 FSSC: http://fermi.gsfc.nasa.gov/ssc/data/ 
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The pulsar was also observed on October 8-9, 2009 with 
the Chandra ACIS-I instrument operated in Timed Exposure 
(Very Faint) mode for a 59.2 ks effective exposure. The Chandra 
Interactive Analysis of Observations software (CIAO v. 4.3; the 
calibration database CALDB v. 4.4.3) was used to generate the 
"level 2" event file for this observation. 



5. Analysis and results 

5.1. Gamma- and X-ray Light Curves 

PSR J1357-6429 is located 2.5° below the Galactic plane 
where the diffuse gamma radiation is still intense. For the tim- 
ing analysis, we took photons within a radius of 0.8° with 
respect to the radio pulsar position (Z = 309.922705°, b = 
-2.514079°) using an energy-dependent cone of radius 6*68 ^ 
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(5.3° x(. 



-)-o.745)2 + Q.09 2 . This choice takes into ac- 



100 MeV> 

count the instrument performance and maximi zes the signal - 
to-noise ratio over a broad energy range dAtwood et al.. 2009). 
The arrival times of events were corrected to the Solar System 
Barycenter usin g the JPL DE405 Solar System ephemeris 
(Standish, 1998), and the events were folded using the radio 
ephemeris from Parkes, using the TEMPO! fermi plugin. 

A total of 8707 gamma rays remain above 100 MeV in our 
energy-dependent circular region, and among them 915 ± 161 
pulsed photons after background subtraction. We estimated the 
background level represented by the dashed horizontal line 
in Figure Q] (top panel) using two circular regions 1.5° from 
PSR J 1357-6429 and at the same Galactic latitude. Nearby 
sources were avoided, and we normalized to the solid angle of 
the source region. An H-test value of 89.6 is obtained above 100 
MeV, corresponding to a pu lsed detection significance of ~8 <x 
(de Jager & Busching, 2010). Figure[T]shows the 25 bin gamma- 
ray light curve (top panel) as well as the 1 .4 GHz radio profile 
(bottom p anel) derived from obser vations with the Parkes radio 
telescope (Wel tevrede et al.Ll2010l) . The gamma-ray peak is off- 
set from the radio pulse by 0.37 ± 0.03 ± 0.01 according to a 
Lorentzian fit with a full width at half maximum (FWHM) of 
~ 0.25. The first phase uncertainty arises from the gamma-ray 
fit. The uncertainty in DM of 0.5 pc cirT 3 causes an uncertainty 
in the extrapolation of the radio pulse time-of-arrival to infinite 
frequency of 0.01 in phase. There is no significant evolution of 
the gamma-ray peak with energy as can be seen from Figure Q] 
(middle panels). We defined the 'off-pulse interval' as the pulse 
minimum between 0.8 and 1 . 1 in phase. 

Figure Q] also shows the pulse profile obtained from the 
XMM-Newton EPIC-pn data with 6 ms time resolution. For the 
timing analysis, 1543 events were extracted (with 270 back- 
ground counts) from a 15"-radius circle centered at the pulsar 
position in the 0.5-2 keV range to maximize the signal-to-noise 
ratio. The photons plugirQ TEMP02 package was used to assign 
a phase to each selected photon. This allowed us to use the full 
complexity of the ephemeris parameters including those for the 
glitch mentioned in Section 3 in the timing model. Performing 
an unbinned maximum likelihood fit of the phases to a constant 
plus single Gaussian we obtained an intrinsic pulsed fraction of 
0.56 + 0.13 (accounting for the contribution of the diffuse emis- 
sion surrounding the pulsar; see Section 5.2), with a Gaussian 
of width 0.48 ± 0.18 centered at phase 0.83 ± 0.03. The H-test 
yields a value of H — 52.42 corresponding to the probability 
of a chance occurrence in a single trial p = 1.6 x 10" n , or 




2 http : //www . physics .mcgill . ca/~aarchiba/ 
photons_plug . html 



Fig. 1. Top panel: Phase-aligned histogram of 
PSR J 1357 -6429 above 0.1 GeV and within an energy- 
dependent circular region as defined in section 15.11 Two 
rotations are plotted with 25 bins per period. The dashed 
line shows the background level estimated using two nearby 
circular regions as described in Section 15.11 Three follow- 
ing panels: Phase histograms for PSR J 1357 -6429 in the 
three indicated energy ranges, each with 25 bins per pulse 
period. Second panel from bottom: X-ray pulse profile ex- 
tracted from the XMM-Newton data of 2009 in the 0.5 - 2 
keV energy band. Two rotations are plotted with 15 bins per 
period. Bottom panel: Radio pulse profile based on Parkes 
obse rvations at a center freq uency of 1.4 GHz with 256 phase 
bins dWeltevre de et all 1201 Oh . 



a 6.7<x detection significance. The main contribution to the H 
value comes from the fundamental and first harmonics, in accor- 
dance with the shape of the X-ray light curve with a single broad 
pulse. Such a pulse profile can be produced by a thermal radia- 
tion emitted from the neutron star surface (see Section 5.2). This 
result unambi guously confir ms the previous evidence of pulsa- 
tions found bv lZavlinl d2007l) in the Chandra data taken in 2005. 

5.2. Spatial and spectral analysis of the X-ray data 

We used the XMM-Newton data from 2009 to derive the spec- 
trum of PSR J 1357-6429. Photons were extracted from a 20"- 
radius circle centered at the pulsar position in the EPIC-pn and 
MOS data. Background was evaluted from similar regions in 
the vicinity of the pulsar. The estimated numbers of source 
counts are 2266 + 72, 1000 + 40 and 996 + 39 for the EPIC- 
pn, MOS1 and MOS2 data sets (respectively). The instrumen- 
tal responses were generated with the SAS rmfgen and arfgen 
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tools. The purely non-thermal interpretation of the pulsar's X- 
ray emission can be most likely ruled out since the single power- 
law (PL) fit to the obtained spectra results in a best estimate 
of the hydrogen column density, hh < 2 x 10 19 cm 2 , signifi- 
cantly lower than the value suggest ed by the puls ar dispersion 
measure, «h - 4 x 10 21 crrr 2 (see IZavlinl d2007l) for details). 
Adding a thermal component, blackbody (BB) or neutro n star 
magnetized atmosphere model (NSA; see IZavlinl (|2009) for a 
review) yields reasonable values of «h and provides good fit 
quality (with xl = 0.9-1; see Figure [2j. Generally, the results 
derived with the two-component model fi t s to th e new data are 
consistent with those reported in IZavlinl (120071) . For example, 
the NSA-plus-PL fit gives the estimate on the neutron star ef- 
fective temperature = (0.95 ± 0.05) x 10 6 K (assuming the 
standard star mass M = 1.4M Q and radius R — 10 km, and the 
distance d = 2.4 kpc). The non-thermal component is well de- 
scribed by the PL model with a photon index F = 1.43 + 0.14 
and (unabsorbed) flux Fx - 2.0 X 10~ 13 ergs crrT 2 s _I in 0.5- 
10 keV (this range is used for the estimates on all non-thermal 
fluxes below), whereas the hydrogen column density is «h = 
(3.9 ± 0.4) x 10 21 cm" 2 . 

In addition to the results yielded by the new XMM-Newton 
observation, the Chandra ACIS-I data of 2009 reveal dif- 
fuse emission both surrounding PSR J 1357-6429 and extend- 
ing 15"-20" from the pulsar in the north-east direction (see 
Figure|3]l. Despite the small number of counts (55 ± 9) estimated 
in the elongated and bent feature (outlined with the dashed con- 
tour in Figure O, its detection is significant at a 6.1 <x level. The 
spectrum of this feature, extracted with the CIAO tool specex- 
tract that yields both the spectra and instrumental responses, is 
fitted with a PL model of F = 1.51 ± 0.28 and unabsorbed flux 
Fx — 0.2 x 10~ 13 ergs crrr 2 s _1 . Thi s confirms the earlier detec- 
tion of a potential PWN reported by IZavlinl (12007b . The diffuse 
emission surrounding the pulsar (denoted with the elliptical con- 
tour in Figure [3] with the estimated number of 224 + 15 source 
counts), is characterized by a PL spectrum of F = 1.48 ± 0.18 
and Fx — 0.7 x 10~ 13 ergs crrr 2 s _1 . A detailed analysis of 
the detected nebul a is outside the scope of this paper but see 
IChang et al.l d201 ll) for more results on the X-ray PWN. 

But it should be noted the diffuse emission is included in 
the extraction region used to derive the X-ray spectra from the 
XMM-Newton EPIC instruments. Therefore, a fraction of the 
non-thermal flux measured from the EPIC data is contributed by 
the PWN. This is clearly seen in Figure [4] showing the pulsar's 
spectrum extracted from the l'.'5-radius circle centered at its ra- 
dio position in the Chandra ACIS-I data (456±21 source counts; 
see Figure [3}. While the thermal component (and the hydrogen 
column density) is well described by the NSA model yielded by 
the XMM-Newton data, the non-thermal component is fitted with 
a PL model of photon index r = 1 .45 + 0.16 and unabsorbed X- 
ray flux Fx 0.9xl0~ 13 ergs crrr 2 s _1 (with^ 2 = 1.1; Figure0). 
We note that the sum of the non-thermal fluxes estimated in the 
Chandra data for the pulsar, the surrounding emission and the 
elongated tail match well that derived from the XMM-Newton 
spectra. 

Another interesting point to mention is that the elongation 
of the detected X-ray PWN, about 60° to North-East, may in- 
dicate the direction of the pulsar's proper motion (see also 
Mig nani et all 1201 ll) . 

5.3. Fermi LAT spectra and phase-averaged flux 

The phase-averaged spectrum of PSR J 1357-6429 was ob- 
tained using gtlike, a maximum likelihood spectral analy- 
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Fig. 2. Spectra of PSR J 1357 -6429 detected with the XMM- 
Newton EPIC instruments (upper panel) fitted with a two- 
component, NSA-plus-PL, model (lower panel). The middle 
panel shows residuals in the spectral fit. 



sis dMattox et all 1 19961) implemented in the Fermi Science 
Tools 1 . This tool fits a source model to the data. The model 
includes the object under study as well as other gamma-ray 
sources in a region extending well beyond the LAT's angu- 
lar resolution. Sources near PSR J 1357-6429 exceeding the 
background with statistical significance larger t han 5 <x are 
extrac ted from the Fermi LAT 2FGL catalog (lAbdo et all 
1201 ll) . Galactic diffuse emission is modeled using the ring- 
hybrid model ringJ2year_P76_v0.fits. The instrumental back- 
ground and the extragalactic radiation are described by a sin- 
gle isotropic component with the spectral shape in the tabulated 
model isotropJ2year_P76jiource_v0.txt. The likelihood tool gt- 
like further requires tables of the Instrument Response Functions 
(IRFs), which are the energy and direction dependent angular 
resolution, effective area, and energy resolution. We used the 
P7_V6 IRFs. These new models and their detailed description 
will be available through the FSSC. 

Assuming a power-law with an exponential cut-off spectral 
shape, the best fit result is obtained for a spectral index of 1 .5 ± 
0.3 ± 0.3 with a cut-off at 0.8±0.3+0.3 GeV and an integral pho- 
ton flux above 100 MeV of (6.5 ± 1.6 ± 2.3) xl0 _s cnT 2 s _1 . The 
integral energy flux is (1.9 ± 0.2 ± 0.8) xl0~ 5 MeV crrr 2 s _1 . 
The first error is statistical, while the second represents our esti- 
mate of systematic effects as discussed below. To check the as- 
sumption of a cut-off energy in the spectrum, we have also fit the 
same dataset with a simple power-law. The spectral model us- 
ing an exponential cut-off is better constrained with a difference 
between the log likelihoods of ~ 4cr, disfavoring the power- 
law hypothesis. Figure [5] shows the spectral energy distribution 
of PSR J 1357 -6429 as seen by Fermi. The Fermi LAT spec- 
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Fig. 3. Chandra ACIS-I image of PSR J1357-6429 (in 0.5-10 
keV) surrounded by an extended structure. The dashed contour 
indicates the jet-like feature elongated to the North-East direc- 
tion. The cross and the 175-radius circle give the pulsar radio 
position. The elliptical contour (with 7" and 4" major and minor 
semi-axes) indicates the diffuse emission surrounding the pulsar 
(minus that located within the circle). 




Energj (keV) 

Fig. 4. Spectrum of PSR J 1 357-6429 detected with the Chandra 
ACIS-I instrument (crosses in the upper panel) and fitted with a 
NSA-plus-PL model (solid curves). The thermal (NSA) compo- 
nent, as well as the hydrogen column density, is the same as in 
Figure 2. The dashed curve indicates the PL component derived 
from the EPIC spectra. The lower panel shows residuals in the 
spectral fit. 



tral points were obtained by dividing the 100 MeV - 100 GeV 
range into 9 logarithmically-spaced energy bins and performing 
a maximum likelihood spectral analysis in each interval, assum- 
ing a power-law shape for the source. These points are overlaid 
with the exponentially cut off power law fit over the total energy 
range. 

Two main systematic uncertainties can affect the LAT flux 
estimation for a point source: uncertainties on the Galactic dif- 
fuse background and on the effective area. The dominant uncer- 



tainty at low energy comes from the Galactic diffuse emission 
since PSR J1357-6429 is located only 2.5° from the Galactic 
plane. By changing the normalizatio n of the Galactic di ffuse 
model artificially by +6% as done in lAbdo et a D (I2010ch . we 
estimate the systematic error on the integrated flux of the pul- 
sar to be 50% below 500 MeV, 20% between 500 MeV and 1 
GeV, and 7% above 1 GeV. The second systematic is estimated 
by using modified IRFs whose effective areas bracket those of 
our nominal IRE These 'biased' IRFs are defined by envelopes 
above and below the nominal energy dependence of the effec- 
tive area by linearly connecting differences of (10%, 5%, 20%) 
at log(E) of (2, 2.75, 4), respectively. We combine the errors in 
quadrature to obtain our best estimate of the total systematic un- 
certainty at each energy, and propagate through to the fit model 
parameters. 

As can be seen from Figure |5j no significant emission is de- 
tected above 4.5 GeV at the position of PSR J1357-6429. In 
very high energy gamma rays, the H.E.S.S. experiment has de- 
tected significant emission close to the pulsar and significantly 
extended with an intrinsic Gaussian width of 0.20° + 0.02°. 
As a first search for unpulsed emission from the TeV source 
HESS J 1356-645, we fitted a Gaussian of 0.2° to the data in 
the energy band 4.5 - 100 GeV where no pulsed emission is de- 
tected. No significant signal could be observed and we derived 
95% Confidence Level (CL) upper limits on the flux in the 4 
logarithmically-spaced energy bins between 4.5 and 100 GeV 
assuming a power-law shape for the source with fixed spectral 
index of 2. In a second step, we used the off-pulse data selecting 
photons in the 0.8 - 1.1 phase interval and fit the same Gaussian 
of 0.2° to data in the whole energy band 100 MeV - 100 GeV. 
Again, no significant emission could be detected and we derived 
95% CL upper limits on the flux in the 9 logarithmically-spaced 
energy bins (assuming the same spectral index). These upper 
limits are presented in Figure [5] 



6. Discussion 



iJohnston & Weisbergl ((2006) measured the polarization profiles 
for 14 young pulsars at Parkes. For PSR J 1357-6429 the po- 
larization position angle *P increases with a constant slope of 
^ = 2.7 across the radio pulse. In the simplest RVM sce- 
nario (Rotating Vector Model, Radhakrishnan & CookeL 11969). 
sin a = ^ sin /?, hence the angle j3 between the magnetic axis 
and the line-of-sight is small, f3 < 22°. The lack of an inflection 
in ¥ vs <p and the narrow <p range of the pulse make a full RVM 
fit unconstraining. 

PSR J1357-6429's gamma-ray profile resem- 
bles that of PSR J2229+6114 (lAbdo et all |2009e), 
PSR J0248 +5832 dTheureau 
J1718-3825 dWeltevrede et al 



et all 
2010) 



2011) 



and PSR 

j, IZ.U1UD in having a single, 
wide peak offset from the radio p ulse by a bit less than a 
half -rotation. The "Atlas" update by iRomani & Wattersl d2010l) 
provides gamma-ray profiles calculated for different emission 
models, over a coarse grid of inclination angles a and £ = a + /?, 
and for 4 gap widths w approximated by the efficiency, 
T] = Ly/E. We find L y = (2.1 ± 0.3 ± 0.8) x 10 34 /n(c//2.4kpc) 2 
erg s- 1 and 77 = 0.01/ n (<i/2.4kpc) 2 above 100 MeV, with the 
beami ng correction factor fa as defined by IRomani & Wattersl 
(2010). Two profiles therein resemble the data. They are for 
w = 0.01, a = 55°, and f = 35° and 50° (adjacent plots in 
the coarse (a, grid). This is consistent with small /3 = £ - a 
suggested by the RVM. For these geometries, the Atlas predicts 
fn — 1, that is, a small correction to L y . This luminosity is 
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Fig. 5. Spectral energy distribution of PSR J1357-6449 in 
gamma rays. The LAT spectral points (in red dots) are ob- 
tained using the maximum likelihood method gtlike described 
in Section IB31 in 9 logarithmically-spaced energy bins between 
100 MeV and 100 GeV. The statistical errors are shown in red, 
while the black lines take into account both the statistical and 
systematic errors as discussed in Section lBTH The magenta solid 
line presents the result obtained by fitting an exponentially cut 
off power law to the data in the 100 MeV - 100 GeV energy 
range using a maximum likelihood fit. A 95 % C.L. upper limit 
is computed when the statistical significance is lower than 3 <x 
(red dots). Blue filled (open) stars represent the 95 % C.L. up- 
per limits derived assuming a Gaussian of 0.2° at the position of 
HESS J 1356-645 in the whole signal (in the off-pulse). 



within the distance-dominated scatter of the L v oc V-E/10 33 rule 
shown by Abd o et al . (2010a), with E in units of erg s . 

The new X-ray data collected from PSR J 1357-6429 with 
XMM-Newton and Chandra in 2009 strongly suggest that the 
bulk of the pulsar flux at soft X-ray energies (below about 
1.5 keV) is of thermal origin, making this pulsar the second 
youngest neutro n star w i th a d etected thermal component known 
at this time (see lZavlirj d2007l) for details). As the X-ray pulsed 
flux shown in Figure 1 is extracted in the 0.5-2 keV range, it 
mostly comes from thermal radiation. It explains the sine-like 
shape of the light curve with the single broad pulse that is char- 
acteristic of the majority of X-ray pulse profiles observed from 
pulsars with dominating thermal emission. It is also consistent 
with the phase shift between the pulses in the gamma-ray and X- 
ray fluxes — the pulses originate at different sites. The thermal 
component is emitted from the neutron star surface, whereas the 
high-energy radiation is generated far from the surface (e. g., in 
the pulsar magnetosphere). Yet the origin of the thermal compo- 
nent remains unclear. It may be emitted from a small hot area 
(polar caps) of a radius of 1-2 km on the pulsar's surface, 

as suggested by the interpretation involving the BB model, 
or from the entire surface as indicated by the fits with the 
NSA models. To discriminate between these two interpreta- 
tions, one should model the pulsed flux, as has been done for 
X-rays detected from a number of mi l liseco nd pulsars (e.g., 
Zavli n & PavlovL Il998t iBogdano v et all [2008). This modeling 
should take into account the effects of the neutron star geome- 
try (orientation of the pulsar spin and magnetic axes with respect 



the line of view), surface magnetic and temperature distributions, 
intrinsic properties (anisotropy) of the thermal emission, and the 
gravitational bending of photon trajectories near the star surface. 
If future X-ray observations provide much better photon statis- 
tics, such an analysis can be supported by phase-resolved spec- 
troscopy and energy-resolved pulse modeling. At present, we 
can only speculate that the detected thermal flux is intrinsically 
anisotropic, as predicted by the atmosphere models, otherwise 
the effect of the gravitational bending would strongly suppress 
the pulsations. 

PSR J 1357-6429 also lies within the extended TeV 

source HESS J1356-6 45, proposed as its associated 

PWN ( Abramows ki et all 1201 ll) . No significant signal is 
detected using Fermi LAT data above 4.5 GeV, nor in the 
off-pulse interval b elow 4.5 GeV. In the cont ext of the leptonic 
model proposed by lAbramowski et al] d201 ll) . our upper limits 
constrain the spectrum of the electron population responsible 
for the extended TeV emission. The best model reproducing 
both the radio, X-ray and TeV data as well as the non-detection 
by Fermi is obtained for a magnetic field of 4 fiG (similar to 
what is seen in other evolved TeV PWNe) and an electron 
spectrum characterized with a spectral index of 2.4, a mini- 
mum energy of 10 GeV, a cu t-off energy of 100 TeV and a 
total energy of 4 X 10 47 erg (lAbramowsk i et all 1201 ll) . This 
energy is only ~10% of the total kinetic energy since birth of 
PSR J 1357-6429, assuming an initial period of 80 ms. 
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